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The exact solution for the thermodynamic and dynamic properties of 
the infinite-dimensional multi-component Falicov-Kimball model for arbitrary 
concentration of d- and f-electrons is presented. The emphasis is on a descrip- 
tive derivation of important physical quantities by the equation of motion 
technique. We provide a thorough discussion of the f-electron Green func- 
*^j ■ tion and of the susceptibility to spontaneous hybridization. The solutions are 

used to illustrate different physical systems ranging from the high-temperature 
phase of the YbInCu4 family of materials to an examination of classical in- 
termediate valence systems that can develop a spontaneous hybridization at 

r = o. 
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Introduction 

The anomalous features observed in the YbInCu4 family of intermetallic compounds 
(Sarrao et al 1999) seem to be driven by the short-range Coulomb interaction be- 
tween highly-degenerate Yb f-holes and the conduction states (Freericks and Zlatic 
1998). The same interaction seems to be responsible for the optical anomalies 
in SmB 6 and related compounds (such as correlated ferroelectrics) ( Wachter and 



Travaglini 1985, Guntherodt et al 1982, Portengen et al 1996). We study the ef- 
fect of this interaction using the multi-component Falicov-Kimball model (Falicov 
and Kimball 1969) in infinite dimensions, where all the thermodynamic properties 
can be calculated exactly. The model consists of (2s + l)-fold degenerate mobile 
d-electrons and static (2S + l)-fold degenerate f-electrons, which interact via an 
onsite Coulomb interaction U. The model was originally introduced to describe 
metal-insulator transitions in materials that do not change the character of their 
electronic states, but do change their thermodynamic occupations as functions of 
the external parameters. The exact results for the static and dynamic correlation 
functions of the spin-one-half model explain the collapse of the low-temperature 
metallic phase of YbInCu4-like systems, and account in a qualitative way for most 
of their features in the paramagnetic, semiconducting, high-temperature phase. The 
exact solution of the spinless model shows that the statistical fluctuations give rise 
to a logarithmic divergence (in T) of the spontaneous hybridization correlation func- 
tion at zero temperature, so that any amount of quantum mixing could lead to a 
phase transition at finite temperatures. This might be relevant for SrnBg and for 
correlated ferroelectrics (Portengen et al 1996). 

In what follows, we describe the model, explain the method of calculating 
the Green's functions for d- and f-electrons, and present results for some static 
and dynamic correlation functions of a spinless and spin-degenerate case. Detailed 
comparison with the experimental data will be given elsewhere. 

Formalism for the d-electron Green's function 

The multi-component Falicov-Kimball model (Falicov and Kimball 1969) describes 
the dynamics of two types of electrons: the conduction electrons (created or de- 
stroyed at site i by d ig or d ia ) and localized electrons (created or destroyed at site 
% by fj or f irj ). The (2s + l)-fold degenerate d-states and the (2S + l)-fold de- 
generate f-states are labeled by a and rj, respectively. The multi-component model 
is used to describe the electrons with spin and/or orbital degrees of freedom, and 
2s + 1 and 2S + 1 can assume different values. The non-interacting conduction 
electrons can hop between nearest-neighbor sites on the D-dimensional lattice, with 
a hopping matrix —tij = —t*/2yD] we choose a scaling of the hopping matrix that 
yields a nontrivial limit in infinite-dimensions (Metzner and Vollhardt 1989). The 
f-electrons have a site energy Ef, and a chemical potential \x is employed to conserve 
the total number of electrons nd + rif = n to t- The d- and f- number operators at each 
site are n d = J2 a n do . and nf = Ylri n fv We assume an infinite Coulomb repulsion 
between f-electrons with different label rj, and restrict the f-occupancy at a given 
site to ri/ < 1, regardless of the degeneracy of the f-state. The Coulomb interaction 
U between the d- and f-electrons that occupy the same lattice site is finite, and the 
Falicov-Kimball (FK) Hamiltonian for the lattice is defined as (Brandt and Mielsch 



1989, Freericks and Zlatic 1998), 

H F k = ^2{-Uj - H&i^Aedjo + ^2(E f - n)fj v f iv + U^ 4*divf! n fiv (!) 

ij,cr i,r] i,cv 

The FK lattice model ([!]) can be solved in infinite dimensions using the methods 
of Brandt and Mielsch (Brandt and Mielsch 1989). We consider two kinds of lat- 
tices: (i) the hypercubic lattice with a Gaussian noninteracting density of states 
p(e) = exp[— e 2 /t* 2 ]/(y / 7rt*); and (ii) the infinite-coordination Bethe lattice with a 
semicircular noninteracting density of states p(e) = ^/4t* 2 — e 2 /(2nt* 2 ); and we take 
t* as the unit of energy (t* = 1). The method of calculation is formulated for arbi- 
trary values of s and S labels, but the results are presented for the spinless model 
(2s + 1 = 1 and 2S + 1 = 1), and for the spin-one-half model (2s + 1 = 2 and 
2S + 1 = 2). We consider only the homogeneous phase, where all quantities are 
translationally invariant. 

Mapping onto the Falicov-Kimball atom 

Infinite-coordination lattices can be solved by a mean-field-like procedure, because 
the self energy of the conduction electrons is local (Metzner and Vollhardt 1989). 
That is, the Dyson's equation for the local d-electron Green's function Gf oc (z) on 
the lattice reads 

where z is a complex variable, and S CT is the momentum-independent self energy. 
Hence, as noted by Brandt and Mielsch (Brandt and Mielsch 1989), the lattice self 
energy coincides with the self energy of an atomic d-state coupled to an f-state by 
the same Coulomb interaction as on the lattice, and perturbed by an external time- 
dependent field, A ct (t, t'). For an appropriate choice of the A-field, the functional 
dependence of E°" on G ff (z) and F CT (z), the atomic propagators for d- and f-states, is 
exactly the same as in the lattice case. The lattice problem is thus reduced to finding 
the atomic self-energy functional for d-electrons, and then setting Gf oc (z) = G a (z) 
and Fi oc {z) = F a (z) at each lattice site. 

The FK atom can be solved by using the interaction representation, such that 
the time dependence of operators is defined by the atomic Hamiltonian, 

H at = -n Y, did, + (E f - n) Y flU + UJ2 did, f v U (3) 

a r) ar\ 

and the time- dependence of the state vectors is governed by the evolution operator, 
S{\) = T T e-f° dT f° dT '^ x ° {T ' T ' )dl{T)d ° {T,) . (4) 



The external field A°"(t, r') is assumed to be time-translation invariant in (imaginary) 
time, 

\°(T,r')=TY,Ke~ lU)n iT - T '\ (5) 



and hence can be expanded in a Fourier series in the Fermionic Matsubara frequen- 
cies u n = 7r(2n + 1)T, where we set ks = 1. In the absence of an external magnetic 
field, the A-field is the same for all the a-components. The unperturbed atomic 
Hamiltonian @ conserves the number of f- and d-electrons, while the time depen- 
dent A-field gives rise to fluctuations in the d-occupancy. In the equivalent lattice 
problem, the local d-fluctuations are due to the d-electron hopping. 

The thermodynamics of the FK atom follows from the partition function, 

Z at {\) = Tr d/ [T T e-^*3{X)] , (6) 

where the statistical sum runs over all possible atomic configurations, which is de- 
termined by function A CT (r, r') for t,t' G (0,/?). The specific feature of the atomic 
FK model is that the number of f-electrons is a constant of motion, and is either 
zero or one, while the d-electron number can assume any value between zero and 
2s + 1. Furthermore, the evolution operator S(\) never transfers the state vectors 
out of the invariant (n/ = or rif = 1) Hilbert subspaces, so the matrix elements 
can be evaluated within each invariant subspace by replacing the operator ^2 f\fr) 
by its eigenvalue or 1. The trace in (|5]) can thus be performed separately for the 
rif = and rif = 1 subspaces. Within the rij = subspace, the operator dynamics 
is governed by a simplified non-interacting Hamiltonian, 

flo = -/*X)4<k, (7) 



and we have dj-(r) = dj.exp(— /ir) and d a {j) = <i CT exp(/ir), where d) a = dj.(0) 
and d ff = rf CT (0) are the time- independent Schroedinger operators. The operator 
dynamics in the n/ = 1 subspace is governed by the same Hamiltonian as in the 
rif = subspace but with \x replaced by \x — U. The trace over f-states in (|B|) can 
be performed, and we find, 

Z at (X) = Z (A,/i) + (2S + l)e^ E f-^Z (\,fi- U), (8) 

where 

Z (A) = n^(A CT ), (9) 



and 

Z°{\°) = Tr d [T T e-^S(X')] . (10) 

The factorization @ holds because the time evolution due to H is such that the 
operators with different a-labels commute regardless of their time arguments, and 
the S-matrix (^) does not change the a-label of a given state vector. Thus, the 
Hilbert space can be decomposed into invariant a-subspaces and the trace in (fTOJ ) 
is over the non-degenerate d CT -states. In each of these subspaces, the operator dy- 
namics is defined by Hfi , and the partition function Z^(A CT ) describes the statistics 
of a d CT -electron subject to an arbitrary time- dependent A°-field. 

In the presence of the magnetic field, the Hamiltonians (H) and ([3D have to be 
supplemented with Zeeman terms, 

Hz = 9d^B Y^ ad ° d ° + 9f^s ^ Vflfv ( n ) 



where ga and gf are the g-factors of d- and f-electrons, respectively. The solution of 
the model for H ^ is a straightforward generalization of the H = case, which 
is presented here (the main difference is that the effective chemical potentials now 
have a spin index dependence). 

Generalized partition function for the Falicov-Kimball atom 

The self-energy functional for the FK atom is calculated by using the equations 
of motion (EOM) for the Green's function obtained by functional differentiation of 
the generalized partition function Z at (X) (Kadanoff and Baym 1962). Equations @ 
to fllCf ) show that we can find Z at (\) by solving a statistical problem for a single 
non-degenerate d CT -state coupled to the periodic A^-field. Consider the contribution 
to Zq^X 7 ) due to the shift of the A-field from an initial configuration A CT (r, r') to the 
final configuration A CT (r, r') + S\ u (t, t'), 

5Z° = Tr d (T T e-^5S(\°))- (12) 

The 5S(\ cr ) is obtained from the usual rules of the calculus of variations, 

6S(X a ) = Sexpi- f dr f dT'A ff (r,T / )4(r)d ( r(r / )} 

rP rP 

= -S{\) dr t/r , 5A CT (r,r')4(r)t/ (T (r , ), (13) 

Jo Jo 



and the time-ordering is taken into account when this result is substituted into (|T2|). 
Performing the substitution gives 

rP rP 

5\nZ° = dr dr'8\' T {T,T')G°{T',T), (14) 

Jo Jo 

where 

G° (t,t>) = -^Tr d {T T e-^d a (r)dt(r')S(X)} , (15) 

^o 

is the d-electron Green's function restricted to configurations with no f-electrons. 
The function multiplying the variation 5X (7 (t',t) is, by definition, the functional 
derivative of the partition function, 

G ^ r) = -6x^yy (16) 

Expressing Gq(t, t') in Eq. fli"4]) in terms of its Fourier components, 

GUr,r')=Tj2G^e- lulAT - T ' ] . (17) 

n 

and using (||) for S\ a (r, r') leads to the functional relation, 

6\nZ^ = J2GoJK (18) 

n 

where G^ n is now defined as a simple partial derivative, 

_ _mnzs_ 
G °"~ a% ■ [ ' 

An explicit expression for GQ n (\°) would allow us to obtain Zq by solving the 
differential equation fll9]) . Functions G^ n and A£, or Gq(t,t') and A cr (r',r), can be 
considered as matrix elements of integral operators Gq and A ", and Eqs. ( |I~4] ) and 
fll9|) can be written in the operator form as, 



<5 1nZ CT = Tr {G%6\°} , (20) 

where the trace denotes an integration over time if we are using non-diagonal ma- 
trices in the r-representation, or a Matsubara summation, if we are using diagonal 
matrices in frequency space. 



Our next step is to find an explicit expression for Gq u , and solve ( |i"9|) to find Z£ . 
The Green's function in Eqs. ( |T5| ) and (|TUp are obtained by the EOM. Consider first 
the case r > r' . To compute d/drG^r, r'), we must first compute the derivative of 
[T T <i -(r)(iJ.(r / )S'(A <J )] with respect to r. It is important to note that the differential 
operator does not commute with the time-ordering operator, so one must proceed 
carefully. Note that when we take a derivative with respect to S(X a ) it will bring 
down terms like d^(r)d a (T 2 ) or dl(Ti)d a (r), and the latter terms will not contribute 
when multiplied by d a (T) — that is, the time-ordering with respect to T\ is the only 
important variable to consider when taking the derivative. So we write the full time 
ordered product in two pieces 



T T d (T)4(/)S{X ff ) = [T T S(A ff )]d ff (r)[T T 4(r')5(A ff )], 



(21) 



with 



S(\ a ) = exp 
S{X a ) = exp 



f dn [ dT 2 X (T 1 ,T 2 )4(T 1 )d IT (T2) 

T JO 

f dn [ rfr 2 A CT (r 1 ,r 2 )4(r 1 )d CT (r 2 ) 

o Jo 



(22) 



Now the derivative can be computed directly to yield 



d_ 



T T [d a (r)di(r')S(X^} = T T 



{nd a { 



dT 2 X ff {T,T 2 )d a (T 2 )}dUr')S{X a ) 



(23) 



where we employed identities like (i -(T 2 )(iJ.(r)(i -(r) = d a (T 2 ) for r 2 > r (which can 
be easily derived from the fact that the time dependence of the operators involve 
only an exponential factor and the anticommutator of the Fermionic operators is 
1). Since the form for the Green's function is different for r < r', one must repeat 
the derivation there (with the same result). Hence the Green's function satisfies the 
following EOM, 



{— + ai)GJ(t, t') - J P dr"X°(T, t")G° (t", r') = 8(t - t>), 



(24) 



where the (^-function arises from the discontinuity in Gq at r = r' . This EOM can 
also be written as, 



li 



dr"[G^]- 1 (r,r")GUr",r') = 6(T-r'), 



(25) 



or, in matrix representation, 

\GlVGl = 1, (26) 

where [Gq] -1 is the non-diagonal integral operator defined by its matrix elements, 

[Go]-\r, t>) = (-A + h)8{t - t') - A CT (r, r'), (27) 

and the unit operator 1 has the matrix elements 5(r — r'). The Fourier transform 
of Eq. (^7]) gives the matrix elements of [Gq] _1 as, 

[G° ]- 1 (iiu n )=iLu n + ii-\l (28) 

and of its inverse [Gq] as, 

G« Q {iu n ) = — — -. (29) 

lL) n + fl - \° 

The diagonal form of [Gq] _1 and [Gq] in Fourier space is the consequence of the 
time-translation invariance of A ct (t, r'). Using Eqs. (0) and (p9|) we obtain the 
differential equation (Kadanoff and Baym 1962), 

1 - 8 ^ (30) 



iu n + fi - \° d\° 

which has to be solved together with the initial (A = 0) boundary condition, 

ZZ(0,fi) = l + e^. (31) 

The partition function for the simplified atomic problem is thus obtained as (Brandt 
and Mielsch 1989), 



Zl (A CT , /i) = 2e^ 2 FT ^ n + /i K . (32) 

J- J- it , I 



and the complete partition function for the FK atom can be written as, 

Z at {\) = \[ZZ{\°,n) + (25 + l)e-K E '-til[Z°(\*,n- U), (33) 



Dynamics of the atomic d-state 

The renormalized d-electron propagator can be obtained, in complete analogy with 
Eqs. (p~2|) - (|16D , as a functional derivative of Z at with respect to the external field, 



such that 



G: t (r, r>) = -I Tr df {T r e-^d a {r)dl{r')S{X^)) . (35) 

Zat 

The difference with respect to the rif = case is that the trace extends now over 
the d- and f-states, including the a and r\ labels, and the statistical operator is the 
full atomic Hamiltonian H at rather than Hq. On the imaginary frequency axis we 
still have, 

G « = --^f' (36) 

which gives, using Eqs. (H) and (0), the result 

where w = Z /Z at and W\ = 1 — w . The weight W\ gives the f-occupation number 
(Brandt and Mielsch 1989). 

On the other hand, starting from the definition of the Green's function in 



(35), and making the usual Feynman-Dyson expansion with U as the expansion 
parameter, we obtain the standard Feynman diagrams, in terms of the unperturbed 
propagators Gq h (t, t'). The self-energy function of the FK atom S^ on the imaginary 
frequency axis is defined by the Dyson equation, 

K = [GlrT 1 - [GT 1 - (38) 

Eliminating Gq u , and hence the A-field, from Eqs. ([37]) and (|38]) yields, 
_ Wl U 



which can also be written in the form given by Brandt and Mielsch (Brandt and 
Mielsch 1989), 



K=l\U-^-±\(^--u) + *u,x£-\. (10) 

* I (Jnn V \ Un/r / VJTt: 



no- / v -' no- 



where the sign of the square root is chosen to maintain the proper analyticity prop- 
erties of E. 

To clarify the physical meaning of the self-energy (|39|) we now perform the 
calculations directly for the original lattice model ([!]) using the EOM. Since we also 
have to consider the Green's function of higher order in the creation and annihilation 
operators, it is convenient to introduce the compact Zubarev notation, where the 
Fourier transformed quantities are written as, 

«A5»„= / d(r~r')e^ T - T '\T T A(r)B(r')), (41) 

Jo 

and (. . .) denotes the thermal averaging over all the states on the lattice, 

(0) = iTr^e-^O. (42) 

The Fourier transform of the lattice Green's function at site % is now written as 
Gna =<C d ia , d\ a ^>„, and the EOM reads, 

(ji + iu n ) < d ia , 4r >n= ! - J2 U i ^ ^' CT ' d * >n +U J2 < dirf infill 4x >« -( 43 ) 



Using the EOM for the higher order Green's function on the r.h.s. of Eq. fl4"5|), and 
considering also the time derivative with respect to the second time-variable r', we 
obtain 

(/i + iu n ) < d ia fl v f iVJ d\ a > n = (fljfir,) - Y^ *v < d ^fi V fv^ d )a >n 

3 

+ U<^d la flf iv ,flf iv dl^> n . (44) 

Here, we used the fact that a given site i cannot be occupied by more than one f- 
electron, which means that the same spin index rj appears for the f-electron operators 
at both times r and r'. As the f-occupation per site is conserved, we also have the 



10 



relation 

< d i(T flfi V , flfivdl >«=< d i(T flf iv , 4 >„ . (45) 

Defining the local (site-diagonal) self-energy E^ as, 

K < d i(T , d) a > n = C/ J] < d itT flf iv: d) a >„ (46) 

we obtain from (|44]) , 

S^ I (/x + io; n ) < d ilT , d\ a >„ + ^ Uj < c^, <ij CT >n ] = 

= UY,(flfiv) + UK « d la , dl »„ (47) 

Denoting the total average f-occupation per site i by w\ = ^2 v (fi V firi), aud using 
(/i + z^ n ) < d iff , d] a > n + ^ Uj < d*, 4 CT >n= 1 + K < 4r, 4 >, (48) 



which follows from Eq. (f43|), we find, 

K (l~(U- ££) « 4r, dl » n ) = C/ Wl , (49) 



which is equivalent to Eq. ( ]5§D and which we recognize as the standard Hubbard-III 
(CPA) self-consistency equation for the self-energy The fact that the exact d- 
electron self energy is given by the CPA, which becomes exact in the limit of infinite 
dimensions for disordered systems, has a simple physical interpretation: As the f- 
electron number per site is conserved, the d-electrons "see" an effective disordered 
alloy potential because at a certain site i there is either the on-site potential 0, if 
the site is not occupied by an f-electron, i.e. with probability w , or there is the 
on-site potential U, if the site is occupied by an f-electron, i.e. with probability 
Wi = 1 — Wq. However, the self-energy functional depends explicitly on u> 1; which is 
not known unless one calculates the full partition function. 
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Formalism for the f-electron Green's function 

The atomic f-propagator, F(t), is defined in the interaction representation for r > 
as (Brandt and Urbanek 1992), 

F v (r) = -Tr d/ {e-^-5(A)/,(r)/t(0)} /Z at (X), (50) 

where the trace is taken over the atomic d- and f-states, including the spin labels, 
and H at is the full atomic Hamiltonian defined by Eq. (§). The time evolution of 
the f- and d- operators is now defined as, /^(r) = e rHat f T1 e~ rHat , and d\{r) = 
e rBat d\ e~ rHat , which leads to the EOM, 

^4(r) = l-V + Uj2ft(T)fv(r)}dl(T), (51) 



±f„{T) = [-(Ef-fi) - Uj2dUr)dAr)]f v (r). (52) 



The integral representation of ([H]) and (|52|) is, 

4( r ) = e -Mr Tre t/E,/o T ^"4(-")/,(-")^ 5 (53) 

/„( T ) = e -(Ef-rir Tre -UZ,,Kar u 4(r")d„(r")f^ (54) 

which can be written as, 

/„(t) = e-^-^S'(r)f v . (55) 

Here 5"(r) is the time-ordered exponential, 



S'(r) = T T e W Lj2j o " I 



' ^dT'dT"xr(r',T")dUr')d a (T")\, (56) 



and 

X r (r',r") = t/6(r-r')5(r' -r"), (57) 



/ // . 
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with Q(x) the unit step function Q(x > 0) = 1 and 9(x < 0) = 0. Eqs. (|53|) and 
lead to the expression, 



F v (t) = -e-P'-^TrTn, {e-^S(X)S\ X r)f v fl} /Z at {\), (58) 

which shows that the operator f v fl(0) arid the constraint rij = Y^nfvfn — ^ enm - 
inates all of the n/ = 1 states from the trace in fl58|) . In addition, all of the inter- 
mediate states of the system, obtained by applying S S' to an initial state in the 
nf = subspace, remain in the same subspace. This is because the f-operators in S 
and S' always appear as equal-time pairs, / r |(T / )/ 7? (r / ) and just count the f-electrons 
at time r'. Thus, the statistical problem for F v (r) is restricted to a constant nj 
subspace and we can replace the operator expression ^2 = f^f-q in H at and in S S' 
by its eigenvalue or 1. We use nf = for r > propagation and rij = 1 for r < 
propagation. The f-electron Green's function for r > becomes, 

F,(t) = -e-to-^Tn {e-^S(X)S'(r)} /Z at {\). (59) 

with the statistical operator defined once again by H in Eq. (0) and the statistical 
averaging performed with respect to all possible states of a d-electron perturbed by 
the A-field and by the additional time-dependent field Xr 

In the interaction representation defined by Hq, the time dependence of annihi- 
lation and creation operators is again d a {r) = exp(/ir)d CT and d) a {r) = exp(— [iT)d\, 
and the time-ordering becomes trivial. Thus, the product of two T-ordered expo- 
nentials in the expression (p3), can be written as a single T-ordered exponential, 



SCK) = T T expl-J2J Q 1 






drdr A CT (r,r )S.{r)d a {r )) . (60) 



where 

A CT (r',r") = A ff (r',r")-X T (r',r"), (61) 

which also depends on the external time r. The f-electron Green's function becomes, 

F,(t) = -e-^f-^T T TT d {e-^S(A)} /Z at (\), (62) 

and the problem is reduced to the evaluation of the statistical sum of a single atomic 
d-level coupled to a time- dependent A-field. 
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Effective partition function for the f-electron problem 

The effective partition function required for the f-propagator can be written in the 
factorized form, 

2o(A) = T T Ti d {er^°S(A)} = \{Z°(\ a ). (63) 



The factorization ( |63| ) holds because the time evolution due to H is such that the 
annihilation and creation operators with different spin-labels commute, regardless 
of their time arguments, and the exponential operator S(\) can be factorized. 

The time-translation-invariant component of the A-field is determined by map- 
ping the FK atom onto the FK lattice, while the additional x T -component is defined 
by the Coulomb interaction between the f- and d-electrons during the time interval 
(0, r). The presence of this additional time-dependent field can be understood as fol- 
lows. In the FK atom the operator dynamics is defined by H at and the d-occupancy 
of the state vector is time-dependent because of the A-field. The interaction be- 
tween f- and d-electrons gives rise to fluctuations in the f-level position and makes 
the potential energy of the system time-dependent. In the effective d-electron sys- 
tem described by Eq. flB2|), the changes in interaction energy due to the fluctuating 
f-d interaction energy are represented by the Xr-field. In other words, the f-electrons 
with an infinite lifetime, acquire dynamics due to the coupling to the d-electron fluc- 
tuations. In this respect, the FK problem is similar to the X-ray edge problem(Si et 
al., 1992). However, while the X-ray problem is formulated as a single site problem, 
in the FK atom the self-consistently determined A-field keeps track of all other f-sites 
in the lattice. At high temperatures, where the coherent scattering of conduction 
electrons on the f-ions could be neglected, the single-site model (where there is no 
self-consistency to determine the A-field) might be representative of the lattice(Si 
et al., 1992). But at low temperatures, where coherence develops, the single-site 
description and the X-ray analogy might not be appropriate, and the lattice effects 
due to the A-field have to be taken into account. 

The partition function Z (X) cannot be calculated with the same procedure as 
for Z (X) because the A-field is no longer a function of r — r'. Thus, the effective 
propagator associated with Z (X) cannot be diagonalized by a Fourier transforma- 
tion, and a simple differential equation for Zq(X) cannot be derived. Nonetheless, we 
use the functional differentiation of Zq(X) with respect to A to generate an effective 
Green's function, 

«.(^ = -|^l (64) 
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such that, 

9Ar',r") = --^TT d (T T e~^d a (r')dUr'')S(K)). (65) 

Similarly, we introduce an auxiliary Green's function defined for a d-level 
driven by the Xr-field only, 

_ <5rnZ CT ( Xr ) 
^ (r ' r) - - 5x T (r",rO (66) 

where 

^( XT )=r r Tr d {e-^5( Xr )}, (67) 

and 

g 0tT (r', t") = ~^r^U {T T e^ H od ff {r')dl{r")S{x T )) . (68) 

The Green's functions p CT and g§ a depend explicitly on times r' and r", and implicitly 
on r. 

The evaluation of Zq(xt) an d (?oct is straightforward because the evolution 
operator S(\t) does not change the number of d-electrons. The Hilbert space for 
d-states, in the absence of the A-field, comprises only two states (rid = and n^ = 1), 
and the result for the partition function is simply, 

Z°{xr) = l + e^- UT . (69) 

and 

Zo(Xr) = l[ZZ(Xr). (70) 

cr 

The time-ordered product in Eq. (|6g|) has to be treated with some care, because the 
functional form of x r -field is not the same in all the parts of the (r , r )-plane [see 



/ // X 



Eq. (|57|) 1 . Eventually, we obtain the following expressions for go a (r , r ) (Brandt and 



Urbanek 1992), 

^oa + = (eo-l)e ( ^^'e^-^", for r" < r < r, (71) 

9oa- = ^~ U)T ' e"^-^" , for r < r" < r, (72) 
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9ob+ 



(fo - l)e- UT e» T ' e-^- u)T " 



for 



T < T < T 



(73) 



/ // 



9ob- = 4oe e 



Ur Jh-U)t -\jfr 



for r < r < r 



(74) 



/ // 



9oc + = (£o - 1K T e"^ 



// / 



for T < T < T 



(75) 



<?oc- = e e^ T e"^ 



for 



r < r < r 



(76) 



where the symbol ± in the subscript relates to evaluating g (r , r ) above and below 
the line t = t and the letter a, b, or c refers to a different region on the — f3 
square as depicted in Fig. (|1|). Here the symbol £ = 1/(1 + e Ur ~ /3 ^). The function 
9o<t{.t',t") is implicitly time-dependent, because its functional form depends on the 
relative magnitude of r with respect to r' and r" . 



P 



>) 







b- 



b+ 



x' 



P 



Figure 1: Different regions of the integration range for r' and t" . 



To proceed, we consider the periodic A-field as an additional perturbation to 
the statistical problem defined by Z (xt)- Thus, writing the full S-matrix in Eq. (|56| ) 
in the factorized form S(X) = S(X)S(xt), we find that the Green's functions ob- 
tained by the functional derivatives 5 In Z(X)/ 5 X(t",t') and <51n Z (xt) I '^Xt( t " '/) 
are related by a Dyson equation. In the operator form, this can be written as, 



9 a — 9oa — 9o<rX a g a , 



(77) 



or, equivalently, 



g- l = g^-X° = g^(l-g Qa X°) 



(78) 



where g a and go a are non-diagonal integral operators both in the time and in the 
frequency representation, while X a is diagonal in the frequency representation. Next, 
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we show that the partition function can be written as, 

Z '(A') = det|^ 1 |. (79) 

This holds, because the functional derivatives 8\iiZq(\ )/8\ (t 1 , r) and 8 In Zq (X a ) / 8 X a (t' , r) 
define the same Green's function, g(r,r'), so that we can write, 

rP rP rP rP 

8\nZ°{X a ) = - dr dr'g a (T,T r )5X' 7 (r',T) = - dr dr'g^r, r') 8X a (r', r) (80) 
Jo Jo Jo Jo 



where g a (r,T') is given by (pop. In other words, the variation of ln^(A) is not 
changed if A cr (r',r) is shifted with respect to some arbitrary but fixed surface 
[go]^ 1 (r' , t) . Using (ff8|) we write, 8X a = 8(X a — [go] -1 ) = ~8g^ 1 , anc ^ obtain, 



8\nZ°{X) = - ! dr I dr'g^r') 8(X° - g£)(/,T) = [ dr( 9(T 8 g^)(r,r), (81) 
Jo Jo Jo 

where in the last equation we arrived at the diagonal matrix elements of (g^Sg^: 1 ) 
by carrying out the matrix multiplication of g a (r, r') and 8g~ 1 {r' , r). Since g a is the 
inverse of g~ l , Eq. flST]) can be written as, 



8\nZ°(X) = / dr<J Info; 1 ] (r,r), (82) 

Jo 

which shows that 5 \nZ§ follows from the variation of Tr lnfg" 1 ]. Thus, 

lnZ CT (A) = Tr \n\g~\ (83) 

and the matrix identity Tr In A = In det A leads to equation (|79|). The Dyson equa- 
tion, Eq. (^), then yields the effective partition function as a continous determinant, 

Z°(X) = det\g£\ det|l - g 0tT A CT |, (84) 

which can be written as 

ZZ(X„) = ZKxr) det|l - g 0a A CT | = Z CT ( Xr ) detD*. (85) 
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Here, we introduced the integral operator D a , which is defined in the r-representation 
by its matrix elements 



D T 1 ,T 2 = Hri-T 2 )- [ dr"g Qa (Ti, 
Jo 



t»)\°(t»,t 2 ). (86) 



The Fourier transform of (pq) defines the integral operator D a in frequency space 
where its matrix elements form a discrete set. The time-translation invariance of 
the A-field leads to the expression, 

D; q = 5 pq -^- [ dr'e^T' [ dr" g 0a (r' ,r")e~^" . (87) 

P Jo Jo 

where the integrals 

M; g = / dr'e^ T ' f dT"e~^ T "g 0a (T',r") (88) 

Jo Jo 

are given in the Appendix. The implicit dependence of goo-( r , T ) on the external 
time makes Mf_ explicitly r-dependent. Since the determinant of an operator is the 
same in any basis, we can calculate Zq (A) in the Matsubara representation using 



the discrete matrix elements (187]). 

The final form for the f-electron Green's function is thus (Brandt and Urbanek 
1992), 

v{) ~ z at {\) ■ m 

It is useful to examine the expression in Eq. (|8~9p for the limits r — > + and r — > (3~ 
for the spinless case. In the former case we have g (iuj p , iuo q ) = 5 pq /(iuj p + /i) and 

F(T = 0+) = _2e^n„(^ + A-A„)/(^) = _^ (9Q) 

and in the latter case we have g (iuj p , iu q ) = 5 pq /(iu p + ji — U) and 
= g-y. = e-^f-^2e^- u »> rU^n + /x - A n - U)/(iu n ) = _ w 

1 } z at {\) u { ] 

as we must have by directly evaluating the Green's function. 
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Formalism for the spontaneous hybridization 

Recent work by Sham and collaborators (Portengen et al 1996) proposed that the 
Falicov-Kimball model may have a ground state that has a nonzero average for a 
spontaneous hybridization (d/'). Such a state would imply that the Falicov-Kimball 
fixed point is unstable to the periodic Anderson model fixed point. This instabil- 
ity is not allowed at any finite temperature because the conservation of the local 
f-electron number implies the existence of a local gauge symmetry which cannot be 
broken at finite temperature due to Elitzur's theorem(Elitzur, 1975; Subrahmanyam 
and Barma, 1988). This is the same conclusion arrived at from perturbation the- 
ory(Czycholl, 1999) and exact diagonalization(Farkasovsky, 1999). Here we will 
show how to test these ideas exactly in the infinite-dimensional limit. (A mapping 
of the Falicov-Kimball model with a Lorentzian density of states onto the X-ray 
edge problem shows that the hybridization susceptibility can diverge at T = 0(Si 
et al., 1992), but it isn't clear that this behavior will survive when one examines a 
conventional lattice.) 

The static susceptibility for spontaneous hybridization satisfies 

Xhyb = - I dr[Tr df (T T e-^S(X)f(r)d^r)d(0)f(0))/Z at (X) + G(r)F(r)], (92) 
Jo 

for the spinless Falicov-Kimball model. The integrand can be determined by simply 
taking the functional derivative of the f-electron Green's function with respect to 
A(r,0): 



Xhyb = - / dr 
'o 



6F(r) 



+ G(t)F(t) 



13 



.e-( E f-^Z ( X r)6[l-g \] 



^ Zat(X) T S\(T,0) ■ (93) 



8X(r, 0) 
But using the identity det A = exp[Tr In A] tells us 

fe^T =-det[l-g \]{(l-g \)- 1 g } 0tT . (94) 

Substituting this result into the integral then gives 



Xhyb 



[ drF{r) [ dr'[l-g X] 0T l ,g (r',r), (95) 

Jo Jo 



where it should be noted that the auxiliary Green's function g is evaluated with 
the Xt field and therefore must be recalculated for each value of r in the integrand 
(since Xt varies with r). 

Now we introduce a Fourier transform for the r' variable to rearrange this 
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result into the following final form 

Xhyb= drF(T)Ty2[l-g X]^g {iu n ,T), (96) 

ran 

where the Green's function is Fourier transformed with respect to one coordinate 
only 

rP 

g (iu n , t)=T dr'e^'goir', r). (97) 

Jo 

The calculation of Xhyb requires little additional work to what is needed to calculate 
F(t). At each value of r we need only invert the matrix [1 — g Q \] and perform the 
relevant vector product with go and the matrix summation. 

It is interesting to evaluate the susceptibility in the limit U — > 0. Here we 
have F(t) = —e~ T ( E f~^/(l + e~ )3 ( Ef ~^), the auxiliary Green function becomes time- 
translation invariant, so the matrix is diagonal [1 — goX]^ n = (i^m + A t )/(^m + A 4 — 
A m )<5 m „, and the partial Fourier transformed Green function becomes g (iu} n ,r) = 
e UJnT / '(iu n + fi) . The formula for the susceptibility in Eq. (p6|) can now be evaluated 
directly by performing the summation over the Matsubara frequencies to yield 

^=0)=/^) /( '-^f-^ (98) 

where f(x) = 1/(1 + e@ x ) is the Fermi function. As T — > 0, the susceptibility will 
diverge whenever the chemical potential is equal to Ef because the Fermi factors will 
limit the integration to the region e < Ef, which will cause the integral to diverge 
logarithmically. A more careful analysis shows that the susceptibility will behave 
like 

Xhyb{U = 0) -> --p(E f ) InT + constant, (99) 

as T — > 0. Since we expect the susceptibility to be larger in the interacting case, this 
analysis is suggestive that the spontaneous hybridization will continue to diverge for 
nonzero U as well. 

Numerical solutions 

The numerical implementation of the above procedure is depicted schematically in 
Fig. ^ and described as follows: We start with an initial guess for the self energy 
S CT and calculate the local propagator in (g). Using Q58D we calculate the bare 
atomic propagator Gg n and find Z Q and Z at . Next we obtain wq, W\ and find G^ 
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from fl37|). Using G^ n and G°, we compute the atomic self energy and iterate the set 



of equations starting with the new self energy until it converges to the fixed point. 

Self energy £ ^ 
^^ from Eq. (38) ^^v 

Calculate G from Eq. (37) Calculate G from Eq. (2) 

\ I 




Figure 2: Iterative algorithm for determining the self-consistent solution of the 
local Green's function. 

The iterations on the imaginary axis give static properties, like n/, and the 
static spin and charge susceptibilities. Having found the f-electron filling wi at 
each temperature, we iterate Eqs. (fj), (|37|), and fl38|) on the real axis and obtain 
the retarded dynamical properties, like the spectral function, the resistivity, the 
magnetoresistance, and the optical conductivity. At the fixed point, the spectral 
properties of the atom perturbed by A-field coincide with the local spectral properties 
of the lattice. 

The f-electron spectrum and the results for classical intermediate- valence 
materials 

As an example, we consider first the d-electron and the f-electron Green's function 
for the spinless Falicov-Kimball model on a hypercubic lattice, at half filling. The 
interacting density of states for the conduction band is shown in Fig. (|3|), where 
Pd{u) is plotted versus frequency for several values of U in the high-temperature 
homogeneous phase. The interacting density of states is independent of temperature 
at high temperatures but not at low temperatures where the system undergoes a 
phase transition to an AB ordered "chessboard" phase (Freericks and Lemanski, 
2000). A metal-insulator gap opens up at U = 1.5. 

We have illustrated how to calculate the f-electron Green's function, and the 
result for the spinless Falicov-Kimball model on a hypercubic lattice, at half filling is 
plotted in Fig. (f|) as a function of Matsubara frequency, and for several values of U. 
In the limit U — > 0, the Green's function becomes a noninteracting delta-function 
[which behaves like l/(iuj n + fi — Ef) on the imaginary axis]. The sharp rise in the 
Green's function can be clearly seen for small U. We expect that the Green's function 
will become a correlated insulator at the same point that the conduction electron 
Green's function becomes insulating. Unfortunately, the f-electron Green's function 
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Figure 3: Interacting density of states plotted versus uo for several values of U, as 
indicated in the figure. 
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Figure 4: F-electron Green's function is plotted versus Matsubara frequency, for 
several values of U as indicated in the figure. 

is temperature dependent here, and since we are working at finite temperature we 
can only see the gap develop when U becomes large [this is easiest to see by the fact 
that F(iuo) — » 0, which occurs for U larger than about 2.0]. 

Recent calculations on the spin-one-half Falicov-Kimball model on the Bethe 
lattice(Chung and Freericks, 2000) have shown the existence of regions of parameter 
space where the ground state is not a charge-density-wave state or a phase separated 
state, but remains homogeneous down to T = 0. In this region of the phase diagram, 
there are no other competing phases, so the system is eligible to have a T = 
divergence of the spontaneous hybridization. Here we illustrate this behavior for 
the spinless model. Previous calculations have found the possibility of spontaneous 
hybridization to be precluded by other phases(Farkasovsky, 1999; Czycholl, 1999) or 
to occur for "singular" density of states(Si et al., 1992). Here we provide numerical 
evidence for the divergence at T = in a restricted region of parameter space for 
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Figure 5: The average f-electron concentration and the uniform charge susceptibil- 
ity are plotted versus temperature, for several values of U as indicated in the figure. 
Both quantities vary little with the interaction strength. The f-electron concentra- 
tion is the set of curves that decrease as T — > and the susceptibility is the set of 
curves that increase as T — > 0. 

"conventional lattices." . 
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Figure 6: The inverse of the hybridization susceptibility is plotted versus temper- 
ature, for several values of U as indicated in the figure. 

We begin by finding a region of parameter space in the spinless FK model 
where the classical intermediate valence state is stable against phase separation all 
the way down to T = 0. This is simplest to find by repeating the previous analysis 
for the spin-one-half case(Chung and Freericks, 2000) (we perform calculations on 
the Bethe lattice here). We can show that if we choose n tota i = 0.5 and — 1 < Ej < 0, 
then the intermediate-valence state is stable for small U. We choose Ef = —0.75, 
where the intermediate valence state appears to be stable for all values of U (we 
did not perform a Maxwell construction to check for first-order phase transitions). 
A plot of the average f-electron concentration and the uniform charge susceptibility 
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versus temperature is given in Fig. (Q). 

Note how the f-electron concentration remains finite for all U as T — > and 
how the compressibility remains finite as well. The inverse of the hybridization 
susceptibility is shown in Fig. (^). Note how the logarithmic divergence at U — 
is difficult to see in this figure and how the inverse susceptibility decreases as U 
increases. This then suggests that the susceptibility will continue to diverge at 
T = for all U (we know from Elitzur's theorem that it cannot diverge at any finite 
T). 



Results for YbInCu 4 

The numerical results for the spin-one-half FK model exhibit several features that 
one finds in the family of YbInCu4 compounds. We consider here a hypercubic 
lattice with a total electron filling of 1.5 and several values of Ef and U, such that 
Ef > fi(T = 0), since that is the regime where the f-occupation can change rapidly 
as a function of T. 

0.5 




Figure 7: Number of the f-holes plotted versus T for U = 4. The Ef increases from 
top to bottom, and is given by -0.7, -0.6, -0.5, and -0.2, respectively. 

The main results can be summarized in the following way. The occupancy 
of the f-holes at high temperatures is large and there is a huge magnetic degener- 
acy. The f-holes are energetically unfavorable but are maintained because of their 
large magnetic entropy. In Fig. ([?D we show nf as a function of temperature, plot- 
ted for U = 4, and Ef ranging from -0.2 to -0.7. Below a certain temperature, 
which depends on U and Ef, there is a rapid transition from the high-temperature 
phase with a moderate f-occupancy to the low-temperature phase where n/ — » 0. 
The "transition" occurs at a crossover temperature T v and becomes sharper and 
is pushed to lower temperatures as Ef decreases at fixed U. The uniform f-spin 
susceptibility is obtained by calculating the spin-spin correlation function (Brandt 
and Urbanek, 1992; Freericks and Zlatic, 1998) and is given by x{T) = Crif(T)/T, 
where C = g1fi 2 sJ{J + 1)/3/cb is the Curie constant. Thus it is clear that in the 
high-temperature phase, where rif(T) does not change much, the susceptibility ap- 
proaches the Curie law. But as far as the f-d interaction give rise to changes in the 
f-occupancy the susceptibility assumes only the form of an approximate Curie- Weiss 
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law. In the region where nf(T) changes rapidly, the susceptibility exhibits a sharp 
maximum, which separates the magnetic and non-magnetic regions of the phase di- 
agram. Below T v , the f-susceptibility is negligibly small, and the total susceptibility 
is due to conduction electrons and is Pauli like. The other static correlation func- 
tions have also been calculated, and the results obtained in the homogeneous phase 
are discussed in (Freericks and Zlatic 1998). 

The interacting density of states pdioj) for the conduction electrons, calculated 
for U = 4 and Ef = —0.5, is plotted in Fig. (|j) versus frequency for several temper- 
atures. (The zero of energy is measured with respect to p.) The high-temperature 
density of states has a gap of the order of U, and the chemical potential is located 
within the gap. Below the crossover temperature T v , rif is small, the correlation 
effects are reduced, and pd{oj) assumes a nearly non-interacting shape, with large 
Pd(p) and halfwidth W ~ 1. 




Figure 8: Interacting density of states plotted versus ui for U — 4, Ef 
(T v = 0.14), and for various temperatures, as indicated in the figure. 



-0.5 



The intraband optical conductivity is determined by an integral of the spectral 
function (Pruschke, et al., 1995) as 



o~(lu) = u dep(e) / du 



f(u)-f(u + u) 



U) 



A(e,v)A(e,v + u), 



(100) 



where A(e, v) = — -Im [l/(v + p — S(z/) — e)] is the spectral function. The result 
for a{uj) obtained in such a way is plotted in Fig. (||) as a function of frequency, for 
several temperatures. Above T v , we observe a reduced Drude peak around u = 
and a pronounced high-frequency peak around oj ~ U. The shape of u{uo) changes 
completely across T v . Below T v the Drude peak is fully developed and there is no 
high-energy (intraband) structure. However, if the renormalized f-level is close to p, 
the interband d-f transition could lead to an additional mid-infrared peak. 

If we estimate the f-d correlation in YbInCu4 from the 8000 cm" 1 peak in the 
optical conductivity data(Garner et al., 2000), we obtain the experimental value 
U ~ 1 eV. Together with T v = 42 K (Sarrao et al 1999) this gives the ratio U/T v ~ 
200. If we take U = 4 and adjust Ef so as to bring the theoretical value of T v 
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Figure 9: Optical conductivity plotted versus uj/t* for various temperatures. The 
U, Ef, and T v , are the same as in Fig. (|8|). 

in agreement with the thermodynamic and transport data on YbInCu4, we get a 
high-frequency peak in ct(uj) at about 8000 cm' 1 , 6000 cm -1 , and 1500 cm -1 , for 
Ef = —0.75, Ef = —0.7, and Ef = —0.5, respectively. 

From the preceding discussion it is clear that the Falicov-Kimball model cap- 
tures the main features of the experimental data for YbInCu4 and similar com- 
pounds. However, our calculations describe much better the doped Yb systems 
with broad transitions, than those compounds which show a first-order transition. 
The numerical curves can be made sharper (by adjusting the parameters) but they 
only become discontinuous in a narrow parameter range. Our results indicate that 
the temperature- and field-induced anomalies are related to a metal-insulator tran- 
sition, which is caused by a large FK interaction and triggered by the temperature- 
or the field-induced change in the f-occupancy. 

Summary and outlook 

We studied the static and dynamic correlation functions of the infinite-dimensional 
FK model by an equation-of-motion method. The exact solution (Brandt and 
Mielsch, 1989; Brandt and Urbanek, 1992) has been presented for the model with 
an arbitrary number of electrons, and for a (2s + l)-fold degenerate d-state and a 
(2S + l)-fold degenerate d-state. In the large-U limit, and for a range of parameters, 
the spin-one-half model exhibits a transition from a high-temperature semiconduc- 
tor or a semimetal, with well defined f-moments, to a low-temperature Pauli metal. 
The static and dynamic correlation functions show many similarities with the exper- 
imental data on valence fluctuating Yb-compounds, and perhaps on SmB 6 , but the 
crossovers calculated for the spin-one-half model are less sharp than what is seen in 
the experimental data. We believe, the sharpness of the transition in Yb-compounds 
is due to the large entropy change, when Yb ions switch from the magnetic 4/ 13 
configuration (with a 14-fold degenerate f-hole in the J = 7/2 spin-orbit state) to a 
non-magnetic 4/ 14 configuration. The model with a 14-fold degenerate f-state and 
a 2-fold degenerate d-state can easily be solved by the methods explained in this 
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paper, and we expect to see a much sharper transition there. The numerical analysis 
of such a model, and the study of the correlation functions for various values of the 
ratio (2S + l)/(2s + 1), will be the subject of subsequent work. 

A more serious difficulty with the FK model is that it neglects the quantum 
fluctuations of the f-state and considers only statistical fluctuations. That is, the 
lifetime of an f-state is assumed to be infinite, and the width of the f-spectrum 
arises only because the f-electrons couple to density fluctuations in the conduction 
band. Thus, the valence transition in the FK model is accompanied by a substantial 
change in the f-occupancy, and the loss of the magnetic moment is associated with 
the loss of f-holes. In real materials, the loss of the local moment seems to be due to 
quantum fluctuations and lifetime effects, rather than to the disappearance of the 
f-holes, The description of the quantum valence fluctuating ground state would have 
to consider the hybridization between the f- and d-states, and that would require 
a periodic Anderson model with an additional FK interaction. The EOM method 
elaborated in this paper does not produce the solution for such a generalized model. 

The actual situation pertaining to Yb ions in the mixed-valence state might 
be quite complicated, since one must consider an extremely asymmetric limit of the 
Anderson model, in which the ground state is not Kondo-like, there is no Kondo 
resonance, and there is no single universal energy scale which is relevant at all 
temperatures(Krishnamurti et al., 1980). We speculate that the periodic Anderson 
model with a large FK term will exhibit the same behavior as the FK model at high 
temperatures. Indeed, if the conduction band and the f-level are gapped, and the 
width of the f-level is large, then the effect of the hybridization can be accounted 
for by renormalizing the parameters of the FK model. On the other hand, if the 
low-temperature state of the full model is close to the valence-fluctuating fixed point 
with the conduction band and hybridized f-level close to the Fermi level, then the 
likely effect of the FK correlation is to renormalize the parameters of the Anderson 
model. 

This picture is borne out from our examination of the spontaneous hybridiza- 
tion for the spinless FK model on the Bethe lattice. We find that as T is lowered 
the system seems to have a logarithmic divergence in the spontaneous hybridization 
susceptibility at T = 0. Normally we cannot reach such a phase because the system 
will have a phase transition to either a phase separated state or a charge- density 
wave, but we can tune the system so that it remains in a classical intermediate valent 
state down to T = 0. When this occurs, effects of even a small hybridization will 
take the system away from the FK fixed point at low enough T and hybridization 
can no longer be neglected. 
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Appendix 

The matrix elements in Eqs. fl87|) and (|Sq), which are used to calculate the deter- 
minant in Eq. ( p9[ ) for the f-propagator, are given by the expressions, 

M( n ^ m ) = ioA nm + (£ — l)B nm (101) 

M (n>n) = £ C nn + (£„ - l)D nn (102) 

where f = 1/(1 + e^"^), and 
.4 
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i(u m -uj n )(iuj m + fi-U) (iu m + ft) {iuj n + /i) 
(io; m + fi)(iu; n + ll-U) 
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iu;„ + /i - i(cu m - Ci ; n )e^-(^m+M)T 

i(u; m - u; n )(iu; m + //)(ia; n + //) 
g-iK-^jr^ + ^ _ [/ + ^ Um _ Un y{iu m +n-u)r^ 
i(u m - u n )(iuj m + [i. - U){iu n + ll-U) 



{iu n + ll — U)t — 1 + e 
(iu n +fi-Uf 



-(Wn+fJ,—U)T 



(103) 



(104) 
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+ 



(iu n + n)(iuj n + /J.-U) 

(ico n + /i)(/3 - r) - 1 - e -^+(^n+M)r 



(iuj n + /i) 2 

(105) 



D, 



+ 



(iu n + fi-U) 2 
(iu n + jj)(iu n + /I - U) 



+ (^ ra + /i)(r-/J)-l-e^-(-" + ^ 
{iuj n + //) 2 

These expressions generalize the matrix elements obtained in (Brandt and 
Urbanek 1992) for the system with electron-hole symmetry. We checked, that for 
fi = U/2 our expressions agree with those of Brandt and Urbanek (1992) for m ^ n, 
but we find a slight inconsistency for the case m = n. The formulae given here, and 
those given by Brandt and Urbanek agree only if the minus sign in front of the iu m 
term which appears in the numerator of the first line of the expression for M( mm ) 
in Brandt and Urbanek (1992) is replaced by a plus sign. 
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